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SUMMARY 


This report presents the results of a study of solar 
radiation observations in the extreme ultra violet spectrum 

O 

from 200 to 600 A made by the OSO-7 satellite. The results 
of special interest in this work were those made under the 
influence of attenuation by the atmosphere in the 250 to 
500 km. altitude range. 

Using published molecular absorption cross-sections 
o r o 

at 304 A ?T,d 256 A, we were able to validate the Jacchia 
atmospheric model (1971) and showed that a mean exospheric 
of 1050°K was appropriate' for the sunset data. The cross- 
sections obtained by Hinteregger and Hall (1969) were 

O 

shown to be most appropriate at 256 A, whereas those of 

O 

Knight, Uribe and Woodgate (1973) were employed for 304 A, 
Employing the validated atmospheric model we were then able 
to derive values of 6.3 + 0.6 Barns and 6.9 + 3.0 flwrws 

— O 

for the cross-sections of 0 and N 2 , respectively at 285 A 
and 7.6 + 0.8 B a- w tfc and 4.8 + 3.9 Ikw«*s for the respective 

O 

cross-sections at 355 A. 
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SECTION 1.0 
INTRODUCTION 


Wolf Research and Development Corporation is pleased 
to submit this final report on project NAS W-2595 to NASA 
Headquarters. The report describes the final results ob- 
tained in analyzing OSO-7 data to evaluate atmospheric 
absorption cross-sections in 300-500 km. altitude range and 

O 

the 200-400 A wavelength range. 

All the analysis and results pertinent to the con- 
clusions are presented in this report. In Section 2.0 the 
basic intensity analysis is given, followed by the geo- 
metrical and orbital analysis in Section 3.0. Section 4.0 
gives a discussion of the results analysis and conclusions 
are presented in Section 5.0. 


SECTION 2.0 
ANALYSIS DESCRIPTION 


The OSO-7 experimental tapes contain information 
on the date of acquisition of solar intensity data, orbital 
characteristic information and the observed intensities at 
one or two lines in the spectrum. The experimental data 
files are organized into pages or frames, there being 96 
intensities collected over 15.36 seconds for each line 
being studied presented within each frame. The fundamental 
sampling interval was 160 milliseconds and the data appeared 
as photon counts. 


For the purposes of our analysis the section of the 
quantity of data to average before applying an inversion 
analysis is indicated by the requirement for a suitable 
compromise between minimizing random errors and systematic 
errors. The greater the averaging time, the smaller the 
random error, but the larger the systematic error intro- 
duced by non-linearities in the model. In the work re- 
ported here we assumed that the average intensity at the 
midtime of a 5 second data span is sufficient. Thus we 
have averaged the intensity for three sections, each 32 
points long, within each of the frames and printed out as 
shown in Figure 2.1. In the Figure 2.1, the various 
columns are as follows: 

Frame: Experimental Frame Number 

Midtime: Middle time of the experimental frame (seconds 

from start of day. Universal Time) . 

AVEUV1: Average intensity in the medium wavelength 

of the frame. 


AVEUV2: Average intensity in the long wavelength of 

the frame. 


2 


P Q 

*3 S3 

of 

<o 

CJ 

£ > 
s ta 




= -=— FR4M&-- 

MIOT1ME 

— AVfcUv’l™ 

- AVEUV 

2 = 

I 

31 056.659 

Cel 58 

2.C02 

-0 




3 . 3C2 

— C 

3 

31087.573 

7.760 

6.20G 

-0 


311 C2. 938 

17. 154 

1 O. 7Crt 


S 

31115.297 

27,792 

16.729 

-0 


= 3l t 33. 65o 


-0 

7 

3 It 49*01 6 

45.917 

25.40o 

-0 

0 J1164.3/-9 

54.771 

27.5P3 

-0 

9 

31WJ.J30 

56.667 

30.417 

“33.187 

-0 

1C 


73 


11 

3*210.457 

67.510 

34.458 
3b. 698 

-c 



65.63b - 


13 

31241 .176 

69.302 

-6.C21 

-0 

1* 

21256.529 

— — 67.4 79 — 

36. 033 

-0 

15 

31271.6/6 

65.104 

36. i27 

-0 



6— .479 

3 3 • 781 


17 

31302.617 

31317.977 

61*250 

34*594 

33.573 

-0 




19 

31333.336 

57*240 

55-969—= 

3 1 .64 6 
32. 199 

-c 


213 48. 69 9 


21 

31 Jr.a.059 

54.655 

30.000 

-0 

23 

J I J 1 < » . 1 3 

31394.777 

31410.137 

50. 969 

2 46 

29.573 

-c 





25 

J142S.496 

40.956 

28.354 

-c 

26 * 

31440.059 

* — 49.6 30 * 

27.354 

-0 

27 

31456.219 

49. 1 as 

26.427 

-0 






29 

31406. 9JB 

47. 959 

25.49C 

-0 

* ■- 30“— 

31502.297 

— *».656 ~ 

25.490 

0 

31 

31517.656 

40.940 
• —*47-693 — 

24.930 

23.930 

€ 




33 

3! 5 43* 379 

47. 13S 

46.625 

24.063 

23.844* 

0 




35 

31579.058 

47.771 

22.031 

2t.281 

0 





37 

316C9.616 

47.333 

21.823 

0 

3S 

2 1 625. 1 76 

46.323 

20.740 


39 

31640.539 

45.927 

21 .C94 

0 







31671.25ft 
3 1666. 6 1 7 

44.365 

20.563 

0 





43 

31701.577 

44.31 J 

20.1 04 

0 

4S 

31732.699 

44.002 

18.635 

0 

46 



*1 6. 7 1 9* 


47 

31763.410 

42.727 

16.335 

0 






49 

31754.137 

43-940 

1 7.688 

0 




l n* 333* 


51 

31324*859 

43.646 

10.365 

0 





53 

31955. S7a 

43.917 

18*010 

0 




18. 073 


55 

31606*297 

44.1Q4 

17.250 

0 






57 

31917.016 

44.940 

17.917 
18. 563 

0 





59 

J1V47.738 

44. c 42 

19.438 

0 






61 

3 19 75.457 

43.63d 

18.094 

0 


.2 205* 2 67b ja 0*~ 
, 20966 708B3a 0* 
>1987Z«22»yJ O* = 

, IR7/33161U.J 3* 


1C5559136-;J Oj — . 


— 0.1 1935SC71 


16^274<*12cj^*0* * 


l S645Sl499J=-0*- 
197462 709 71/ 0. 


2193206164a 0« 


2622 1 1 01 7aJ 0* — 


30392lK160-» □» - 


-0.47678J2w3’’D 04 
- 3. 46263^4 36*»D 04 
-J. 4684J46 IbJD 04 
-O.4’143J0->95JD 04 
-0 .49 322 355150 04 
— O .604 33416020 04 
-0.60969472470 04 
-3.514731 94650 04 
-0.6 1 95/409330 04 
-3.524271 501 CD 04 
-O.r>231l00 lOnO 04 
-3.5J31J35J4CD 04" 
-0.63743021 730 04 
-0.54 14 >105650 04 
-O. 54533766590 04 
-0.349119 11 0*0 04 ’ 
-3.5527122/450 04 
-0.5361 1773570 04 
-0.6593^6130^0 04 
-0. 33242 7294nD 04 
-3.365*2140210 04 
-0.3a3153427S0 "04 = 
—3.57071 09 122D 04 
— = 3 • 3 733^5557*^0 04 
—0.5/5 23 015610 04 
-0.5773329/660 04 
— 0 . 37921 2401 80 04 
-3.330 J&t £&>$•'»» = 04 " 
-3.33259Ct45-*D <4 
-3.33 336 f t JC70 94 
3 . 5350 j | I SdbD C* 
0.3861 >4 7o6oD 04 
3.666/7447900 
0.58764937630 

0.3ddl49484oD 
0.38032151640 
3.38667927190 
0.3596625 33 :5D 04’ 
3*5394/122330 04 
0.53313137690’ 04“ 
3.63751956450 04 
0.336/132S610 04 
3 .530 3 j 375660 04 
■3.3353264 2020 04 
•3. 5353000 1960 04 
0.39241600040=04- 
3.3dOa^64 7220 04 
9.6 7913/0 3370 04- 
0.37727174140 04 
J. 37520639990 04- 
3.5/23 71 3 6020 04 
0.37961093140-04- 
0.5^834792930 04 
3. 3n3Jl 577020 C4 
0.362*2*25530 C4 
0.359434197*0 04 
3. 536 1879 u I 70 04 
0.352 J38o45;Cr=04- 
3.34/262967*0 04 
9*3*56/5092/0 04 
0.341 7 j21 65JO 04 


-0.36524616000 


—0.36934699320 


-0.3710 J13682U 
0. 371 1218574D 
—0.371 1952440D 


—0.3609 30906 ID 


-0.36691779210 


04 

04 

04 

04 

04 


-0. J319J53629D 


—C.3*» 537-V70 7 ID 


■0.33790624310 


•C. 32953570090 


-0.29929924930 

-0.29540683930 


—0.26646721 5 5D 
-0.26199976100 


04 

0.4 0126533 38D 

03 

-0.24061223790 

00 

04 

“0.4C309760770 

03 — • 

-0.22387739350 

00 = 

C4 

0.4C492594640 

03 

-0.20703340150 

00 

04 

0 • 4 C57* Si 1 5*? 

— • 

-*T. 19^235^3540 

00 ■*— — 

04 

0.40863 17792D 

03 

-0.17424915890 

00 

04 

— 0.41040042170 

03 “ 

-0.15731285760 

00 

04 

0.4 123216335D 

03 

-0.14033815770 

00 

04 

= i — -0.4 1415350300 

03 -* 

— 0. 123 ?259 6440 

00 — - 

04 

0.416191644SD 

03 

-0.10628563620 

00 

04 

=’ — -* 0.41«05242eiD 

03“" 

-0.892262460 80-01= 

04 

0.4 1990 305370 

03 

-0.72163193600- 

-01 

04 

0.42174121 060 

03 * 

—C. 550 65 964080- 

-01 

C4 

0.42379 5229 ID 

03 

—0.37966700740-01 

04 

0.4 2565973590 

03 = 

-0.20866364400-01 

04 

0.42751 101 OOO 

03 

-0.377*62931 60-02 

04 

- 0.4 2934 722240 

03 “ 0 .13309229950-0 1 — 

04 

0.431 40e/7CoD 

03 

0.30397535580-01 

04 

-•* 0.43326551030 

03 

0.47452706250-01 

04 

0. 4 2510 90 3940 

03 

0.64462193020-01 

04 

0.42693543170 

03 * “ 

* 0.61487281310-01 

04 

0.4 JRi89l674tl 

03 

0-969655 12600-01 

04 -=0.4 4059186370 

03 — 

”0.11589104030 

00 

04 

0.44247196920 

03 

0. t 3 27 726 303 D 

00 

04 

0.4 444709 3 SCO 

03 

" 0.14907609410 

00 

04 

0.44625649530 

03 

0.16585552560 

00 

04 

- 0.44801816170 

03 

0.16258080200 

00 

04 

0.44975272720 

03 

0.19923857530 

00 

04* 

0.45I921137C0 

03 

• 0.21533105490 

00 — 

'04 

0.453679592^0 

03 

0.23234793550 

00 

04 

= 0.4 554 096 6 6bD 

03 

0.24878646290 

00 

04 

0.4571 1211450 

03 

0.265X4274100 

00 

04 

0.4b^2Sj002ttU 

03 = 

0.28141624960 

00 

04 

0.4 €09756*5 CD 

03 

0.297597442/0 

00 

04 

0.4C25696091D 

03 — 

“ 0-31367841510 

0 0 * " 

04 

0.46433 J13120 

03 

0.3296S87594O 

00 

04 

— 0.46643682090 

03 — 

” 0.34553870130 

00 - 

04 

0.46011682360 

03 

0. 3 61 3 03602 7 D 

00 

04 

0.46976690960 

03 0.37695470070 

00 = — 

04 

0.47138555650 

93 

0.39249136040 

00 

04’ 

— —0.4 7344 24b 160 

03 0.40790668660 

00 - 

04 

0.47507302860 

03 

0.42316927960 

00 

04 

~~ 0.47667392270 

03 0.43834079770 

0 0 — ™ = 

04 

0.4702406262D 

03 

0.4533570B04D 

00 

04 0.4E024 15 J3GD 

03 — 0.4 6323978150 

00 — - — 

04 

0.46101857120 

03 

0.432976531 BO 

CO 

U4 

"0.4 8336266 730 

03 

”0 .4975621 1990 

00 — 

04 

0.46*37218770 

03 

0.51 199635040 

00 

04* 

— 0.48680846250 

03 — 

- 0.52628154560 

00 — - 

04 

0.40832460000 

03 

0.54040101060 

00 

04- 

0.4S9S0713940 

03 — = 0.55435760370 

00 — — * 

04 

0.49125447650 

03 

0.S6S1SI07Z7O 

00 

04- 

1 ■*' ' 0.4931 l 76514D 

03 

“ 0.58177726360 

00 "■ 

04 

0.49456768600 

03 

0.59522266640 

00 

04 

— 0.49596364480 

03 — 

“*0.60349061020 

00=’ 

04 

0.49736303940 

03 

0.62157758340 

00 

04- 

~ 0.49914520210 

03 “ — 

' 0.63448692530 

00 

04 

0.5 CO 52 44951 O 

03 

0.64720141900 

00 

04- 

*' ^O.b Cl 86895 1 60 

C3 

“’0.65972795310 

09 " 

04 

0.50317608570 

03 

0.67205680090 

00 

04 

0.5 C4 369 81 170 

03=*“ 

= 0.68419511790 

00 

04 

0.50617345010 

03 

C.6961255089D 

00 


FIGURE 2*1. COMPUTER OUTPUT SHOWING AVERAGED INTENSITY AND ORBITAL DATA 


X 

Y 

Z 


Satellite geocentric coordinates in km. at 
the midtime o£ the experimental frame. 


Altitude: Altitude of the satellite from the Earth's 
surface in km. at the midtime of the ex- 
perimental frame. 

Sun Angle: Cosine of the angle between the lines 
joining the satellite and solar center 
to the geocenter. 

The calculation of the position, altitude and the Sun 
angle of the satellite will be discussed in detail in 
Section 3.0. 

The intensity data at a given wavelength X is meas- 
ured by the satellite instrument. We assumed that when 
the path from the satellite to the Sun did not pass through 
any layer below 480 kms . there was no significant attenua- 
tion. Thus this data was accumulated and averaged to give 
the mean intensity per observation for zero attenuation, 


When the path from the satellite to the Sun fell 
below 480 kms. the mean intensity per observation 1^ gen- 
erally fell below I 0X indicating attenuation. The total 
optical thickness along the path through the atmosphere 
could be estimated as: 



4 


When the Sun-satellite line fell below the 300 km. altitude 
it was found that no significant signal was received. Thus 
the intensity yields information about attenuation in the 
260 to 480 km. altitude range. 

We divided the atmosphere in this altitude range into 
22 layers each 10 kms . thick. This was considered to be the 
greatest degree of detail we could reasonably consider when 
using one file of data (one orbit) only. 

The equation used for analyzing the data was 



-X F ij 



( 2 . 2 ) 


where 


T 


Xj 


is the observation of optical thickness at 
wavelength X. 



i.L 

is the path length in kms. through the i layer 
for the observation (see next section). 


iX 


til 

is the mean absorption cross-section in the i 
layer (km ) . 


is made up of two parts, one for absorption by 0 and the 
other for absorption by N 2 . If C n oi» n Ni^ are t * ie mean den- 
sities of 0 and N 2 , respectively in the i*"* 1 layer and 
and Oj^ are the respective cross-sections at wavelength X 
we have that 


a iX 


= n 0i °0X + n Ni a NX 


(2.3) 


5 


Thus the values of ^ contain information about both the 
0 and N 2 densities, as well as their cross-sections. Clearly, 
we neod to know at least 3 of the 4 unknowns in each layer 
before the other parameter can be determined. Therefore, 
even if the cross-sections are known it is impossible to 
unambiguously derive both the N 2 and 0 densities independ- 
ently. 


The procedure we employed in our analysis was to use 
the standard atmospheres of Jacchia (1971) to determine the 
appropriate range of exospheric temperatures where the cross- 

O O 

sections were known (256 A and 304 A) and then to employ a 
least squares analysis to derive the cross-sections cr^ and 
at the wavelengths where they were not known (285 A and 
355 A) . The details of the results will be considered in 
Section 4.0 of this report. 

In deriving the values of a ^ using equations of the 
form of (2.2), it was necessary to determine the appropriate 
weighting factor to be used tvith each equation. The objec- 
tive was to evaluate the variance, a , in and then divide 
the corresponding equation by a. 


Suppose that Iq^ 
over a time T Q ^ and 1^ 
over a time T-^ . Then 


was derived from a total of N 
was derived from a total of 


counts 

counts 


r ox = N 0A 




(2.4) 
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Then 


Wox) ■ "ox r ol 
v <-h> "h T x 2 


(2.5) 


give the variances in and 1^, respectively, since N 


and are Poisson counts. 


It follows that 


OX 


■ »Sx 


Ii V(I,) « N 


-1 


( 2 . 6 ) 


and 


W ’ W ' '"ox "x> 


(2.7) 


Using (2.1) we then obtain the simple result, 


V(t x ) 



( 2 . 8 ) 


2 

which is the required variance, a . 


A scries of experiments we performed on the inversion 
of overdetermined linear equations indicated that the best 
combination of accuracy and computer time requirement would 
be yielded by the IBM systems subroutines package inversion 
programs LLSQ and DLLSQ. In practical applications, however, 
these solution procedures failed to successfully invert (2.1) 
consistently. The reason for this failure is not clear but 
we believe it relates to the substantial random errors 
occurring in the Tj estimates. The program was carefully 
checked to eliminate any possible dimensioning or other 
errors and the fact that the routine successfully inverted 
some data sets, supports the hypothesis that the routines 
were being employed correctly. 

The lack of consistency in producing successful 
inversions seriously slowed our production analysis activity. 

As a consequence we developed a new procedure which has 
proved to be very stable and reliable, while producing 
sensible answers consistent with those generated by success- 
ful LLSQ runs. The procedure is based on the summation of 
sets of equations of the form (2.1) to generate a right 
upper triangular matrix on the left-hand side. This proved 
to be particularly easy to accomplish since the altitude of 
the grazing incidence increased monotonically during sunrise 

a 

and decreased monotonically during sunset. (This altitude 
establishes the number of leading zero F^ values in (2.1).) 
Thus the procedure required merely the summing of sets of 
adjacent model equations. Non-singularity was assured in 
all practical cases by the stipulation that there should be 
at least one grazing incidence in each layer of the atmospheric 
model. This is discussed more fully in Section 3.1. 

Once the left-hand side matrix is right upper tri- 
angular, the cross-section in the top layer can be solved 
for immediately and successive back substitution employed 


8 


to sections within all other layers. In practice, we have 

found that the solution can bo generated in much less time 

2 

than by LLSQ since the procedure requires N multiplications 
2 

compared to N (N+n) for the standard least squares procedure. 

(N is the number of coefficients and n the number of equations.) 
The results generated by the application of the new routine to 
the data analyzed in the previous report gave cross-section 
profiles very similar to those yielded in LLSQ. 



SECTION 3.0 

ORBITAL INTERPOLATION AND GEOMETRIC ANALYSIS 


3.1 ORBITAL INTERPOLATION 


The data tapes contain orbital data in geocentric 
coordinates at one minute intervals as illustrated in Figure 3.1 
of the first quarterly progress report on this project. The 
problem that we considered was that of obtaining the geocentric 
coordinates at the mid points of the experimental frames using, 
the given data at one minute intervals. 


The first step was to compute the nearest given orbital 
data point to the frame mid point. Clearly this could not be 
more than 30 seconds away. We thus knew the orbital elements 
including velocity components at an epoch t and wished to com- 
pute the elements at a time t + At where At was in the range 
-30 seconds to +30 seconds. Let X^(t) be the ith coordinate 
at time t. Cowell's equations for two body motion is 



(3.1) 


where 


and 


In making the analysis we neglected the time variation of R 
since the OSO-7 orbit is not far from circular (eccentricity 
around 0.0165). Thus we obtained from equation (3.1) 


5 ‘ 3 

y is the gravitational constant = 3.986013x10 km 
sec” 2 


R is the geocentric radius 


c £ i=l V > 


1/2 


10 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Tape NUMCER* CAC614 
FR4Me”NuMB*R= to J *— 

- MEDIUM WAVELENGTH*' 


FILE NUMBER* 62 
YEAR*-- 1972" OAY Or-YEAx* 


' TIMS -QF-OAY= 


“3610C SECOft-OI 


31 *'.V3" ANGSTROM S" 


CONG-*AVE.fcN3Trt= 363.10 ANGSTROMS 


« 

e 

8 _ 

“ 


“SATELLITE POS1TION- 

~s a tele rre- velocity- 
-solar- po sition 


-z 


-0.564S6AG567D -04" 
“o .so* 63S2A2eo~or*' 
_ OiS62Sr,S30 1 30—00“ 


"'0.2J6A2S66aj0 -OA- 

^O..! 19121 B* 30-01” 


*■3 .-2Z113riS330-00“ 


“-o.27e6*26aaoo- o* — - (kms.i— — — 

“-0.I79S33SM8O 01 tXAS./SEC.)' 
O.9S9O211B03O-OT (A.UI 


-TAPE NUMBER*' CA0614 ’ " FILE KJMO£R*- 62 ^ ... 

FRAXE NUMBER— 104 YEAR* 1*72 nAY nc vc — 

— - 1 ,, PAY OF YEA RS nO ___ TIME OF OAY* 36240 SECON OS 

MEDIUM WAVELENGTH* 314i90 ANGSTROMS , nK - . . _ A " ‘ “ 

- ANGSTROMS L ONG «A V Ct .CNGTH* J68.10 ANGSTROMS 

r 


SATELLITE POSITION 


SATELLITE VELOCITY 


-0.54017045290 04 
0 .428049993SO 01 


-9..1/25O14044J 04 
-0...9 1940343*0 01 


SOLAR POSITION 

TAPE MJMQER* CA 062 4 
-FRAME-NUMBER*-- 105- 

-MEOIIJW"WAVELENGTH= 


0 .062568*3870 00 


—0.29474S 18740 04 (KMS.I 

— 0.23693 SOO 940 01 (KMS./SEC.) 


-0.221 140S1830 00 


-0.S5B9749I Q'JD-01 CA.U1 


FILE NUMBER* 62 

-YEARp-1972 DAY’OF-YEAK* - 06— 


-314.90- 


‘SATELLITE- POSITION ■ 
3ATELLITe-VELOClTY- 
SOLAR-POSITION - 


--0. 51323357340 04- 
O'. 40934546230— 01" 


TIME OF - O AY*- -36300 SECON OS—— . — 

;N3r«=— - — 348. TO - ANGSTROMS™— 

—0.30741425170-34 — - -0.3094S62BJ6D 04 (KMS. I 

— 0.56920499810-01 f KMivTSCC. - 


-O .96257203820’ OO - 


-0.321 12984770-00 ■ 


0.95892864230-01- IA.UJ 


“TAPE-NUMBER* CA0614 FTLt- NUMBER*— 62 

FRAME NUMBER* 106 YEAR* 1*72 DAY OF YPAp- - ‘ — 

DAY ^ AM - TaX* OF DAY* 36360 SECON OS 


MED IUM WAVELENGTH* 314.90 ANGSTRCMS LONG -AVElcNSTH* 


368.10 ANGSTROMS 


SATELLITE POSITION -0 .483EB6956ED 04 


—0.34032332810 04 


satellite velocity 

SOLAR POSITION 


-0.32271024480 04 (KMS.I 


0.50333234190 01 


— 0*54336013 7pJ 01 


0.9625—5 1 3760 00 


-0.22111917330 00 


TAPE NUMBER- CA0614 FILE NUMBER* 62 

—FRAME— NUMBER* — 107 Y£«v>^ -i^7 ? 0Ar ~ 

-MEOIUM- RAVE LENGTH* 31 4. VO -ANGSTROM S- 

1 ~ - " — X 


—0.208391 5S910 01 
-0.95888277470-01 


OCMS.7SEC.) 

(A.U ) 


nxe-oF— jay* 


—36420— SECON — 03- 


LONG -wavelength* 

■ “ — V 


—368.TO— ANGST ROMS' 


-S a tell i t e-pos i t i oir 


-z 


-0 .43227*40293-04- 


-3.-37263730710-04' 


— S4TELLITE VELOCITY 
—SOLAR POSIT ION 


-0.-33444873680-0* (KNSVJ- 


0.544 835 I 6220-01* — 
- 0.9625782 3710- 00-— — 


—9.310041 d45IO~0l“ 
— -0.221 10 360 7*0 09 


'-0.18269360210 Ol “(KMS.TSEC.I- 

— 0.95883699730-01 (A.UI* 


FIGURE 3.1. COMPUTER OUTPUT OF ORBITAL DATA 
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and then employed the truncated Taylor series 

N 1 d3jt i 5 

X,(t+At) = X.(t) ♦ xlL, — -rri(t) At J (3.3) 

1 1 } 1 j! dt 3 


The computed elements were employed to evaluate the appro- 
priate altitude and sun angle. In computing the altitude 
we allowed for the earth's flatness,, f, using the relationship 

R e = R q (l-f(l + l.Sf) sin 2 <J> (1 - sin 2 <J>/(l + § f 2 ))) 

(3.4) 

where Rg is the distance from the geocenter to the earth's 
surface at latitude <J> and R ft is the equatorial radius, 
f was taken to be Ranged by a ^out 5 kms. over 

the orbit. 

In Figure 3.2 we present the computed altitude as a func- 
tion of time for three orders of interpolation (2, 3 and 4).. 

(It should be pointed out that the lack of continuity seen 
in the N=2 data was the worst of any data analyzed to date.) 

On the basis of these data we have elected to use fourth order 
near circular orbit interpolation for all future runs. 

The sun angle was computed using geocentric solar 
position coordinates obtained by linear interpolation between 
the one r.inutc values extracted from the data tape. A plot of 
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the geocentric solar position coordinates against time indi- 
cated that the errors arising from this procedure would be 
negligible. 


3.2 GEOMETRIC ANALYSIS 

In order to analyze the optical thickness data we 
needed to compute the distance covered through each atmospheric 
layer of our model. The philosophy employed here was that the 
earth was spherical, the radius being that appropriate for 
the latitude of the satellite at the time of data acquisition. 
This was assumed to give a. reasonable representation of the 
atmospheric layer structure while reducing the computer time 
for the analysis to a minimum. 

Figure 3.3 presents the geometry of the problem in the 
plane of the geocenter, the satellite and the Sun. P repre- 
sents the satellite, the solar angle, of the satellite 
being as shown. On the path to the Sun from P, A is, the point 
of closest approach, distant R sin <f> from the geocenter, where, 
R is the geocentric radius to the satellite. Using 
Pythagorus's theorem the distance covered between the Jth and 
(J+l) the layer beneath the satellite altitude is 

2 fV R J+l 2 " R 2 sin 2 4> - * Vrj 2 - R 2 sin 2 <}> 



whereas the corresponding value for layers above the satellite 
altitude is the same, but with the 2 replaced by unity. Rj 
is the geocentric radius of the Jth atmospheric layer. 
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Clearly, if R sin <J> is less than R E + h m i n w ^ crc r e * s 
the (local) earth radius and h the lowest altitude con* 
sidered and $ is greater than * then no signal can be received, 
and the data point is ignored. 

i 

Suppose that R R :'.s the radius of the lowest layer 
boundary crossed by the path. In this case, the distance 
covered by the path in the Kth layer is just 


Vi 


2 vRir - ;r 


• 2. 

sin <t> 


If P is in the m th layer then the total distance covered in 
this layer in the case shown is 


V 



R 2 sin 2 <J> - 



R 2 sin 2 $ 


- R cos<{> 


The formulae given here were sufficient to evaluate the dis- 
tance parameters of equation 2.2. 
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SECTION 4.0 

RESULTS OF THE ANALYSIS 


The cross-sections per km. were calculated from the 
data in each layer of the atmosphere from 260 km. to 480 km. 

The thickness of each layer was 10 km. The total cross- 
sections from each file for the wavelength of interest 
are plotted. These are shown in Figures 4.1 to 4.8,. 

O 

We have fitted the data for 304 and 256 A models 
with different exospheric temperatures. The three solid 
curves shown in the Figures 4.1 to 4.4 are from these models 
and we see that exospheric temperature 1050° models fits 
well with the data at sunset at 304 and 256 A. The total 
cross-sections obtained from Jacchia models used Knight et al 

o o 

(1973) 256 A and Hinteregger and Hall (1969) 304 A absolute 
cross-sections. We see that the data af lower level layers 
at sunrise is consistent with sunset data; however, at 
higher altitudes the sunrise data has large scatter. This 
scatter possibly can be attributed to some interlocking 
problems and possibly atmospheric instabilities while heat- 
ing. Due to lack of knowledge of the cross-sections as a 
function of temperature and the temperature change itself at 
the different altitudes, we elected at present to look at 
the sunset data only. 

We derived our estimates of nitrogen and oxygen densities 
from the best fitted Jacchia model and presented in Table I 

O O 

for 304 and 256 A. Using these densities in turn at 285 A 

O 

and 355 A , we derived the absolute cross-sections of oxygen 
and nitrogen at these wavelengths as shown in Section 2.0. 

The cross-sections thus derived are given in Table II. These 
cross-sections are new in literature as far as we know. We 

O 

thus have the densities N~ and 0 and cross-sections at 285 A 

O ^ 

and 355 A. 
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FIGURE 4.4. TOTAL CROSS SECTION VS. ALTITUDE AT 256 A° 
FOR SUNRISE 256 A° 
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FIGURE 4.7. TOTAL CROSS SECTION VS. ALTITUDE AT 285 A 
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In Figures 4.9 and 4.10 we present the oxygen and 
molecular nitrogen densities versus altitude for different 
exospheric temperatures of the Jacchia (1971) model that 
was employed in our study. These graphs illustrate that 
a significant change in the density and, therefore, the 
attenuation characteristics of the atmosphere can be 
expected with changing exospheric temperature. 0 becomes 
more predominant with increasing altitude, so information 
about N 2 must be obtained from attenuation through the 
lower layers. 


SECTION 5.0 
DISCUSSION 


We have in this project arrived at the goal which 
is to derive the N 2 > 0 densities as a function of altitude 
from 260 km. to 500 km. above the surface of the Earth. 
However, the data obtained at sunrise is very erratic and 
showed no trend whatsoever thus giving us only sunset data 
to deal with. However, the sunset data was adequate to 

o o 

derive the absolute cross-sections at 285 A and 355 A of 
0 and N 2 which is an unexpected by-product. 

There are a number' of possible explanations for the 
erratic data at sunrise. The sunrise data span was much 
longer in duration than the sunset portion and therefore 
promises much more information about the atmosphere. It 
is likely that our description of the atmosphere as varying 
only in the vertical direction is not valid and that varia- 
tions in exospheric temperature around the globe become 
important. There is ample data to justify an extensive 
survey of this problem but we were unable to study this 
problem in depth within the time and funding constraints 
of the project. It is highly desirable to look at attenuation 
data in this region from other spacecraft. This combining of 
all the data available from many spacecraft will lead to a 
high accuracy density profiles and atmospheric model valida- 
tion. 


We thank NASA Headquarters for giving us an opportunity 
to participate in this program. We, in consultation with the 
Principal Investigator and the Technical Officer, propose to 
publish this work in a scientific journal at a later date. 
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Table I. Densities of Molecular Nitrogen and Atomic 
Oxygen From the Best Fitted Jacchia Model 


NO. ALTITUDE DENSITY/CM 3 




MOL. NITROGEN 

OXYGEN 

1 

260 

3. 591xl0 8 

1.676X10 9 

2 

300 

1.048x10® 

8.204x10® 


3 

340 

3.2 

xlO 7 

4.153x10® 

1 

!' 

p 

i: 

4 

380 

1.02 

xlO 7 

2.146x10® 

1 

* 

( 

i 

5 

420 

3.28 

xlO 6 

1.125x10 s 

: 

I 

! 

6 

460 

1.08 

xlO 6 

5.95 xlO 7 

ii 

li 

t 

1! 

7 

500 

3.61 

XlO 5 

3.18 xlO 7 

II 
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Table II. Cross- 


285 

CONSTITUENT 

MOL. NITROGEN 
OXYGEN 
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